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Abstract 

We consider the Hamiltonian system consisting of a Klein-Gordon vector field and 
a particle in IR 3 . The initial date of the system is a random function with a finite 
mean density of energy which also satisfies a Rosenblatt- or Ibragimov-type mixing 
condition. Moreover, initial correlation functions are translation- invariant. We study 
the distribution m of the solution at time t £ IR . The main result is the convergence 
of fit to a Gaussian measure as t — ► oo , where /ioo is translation- invariant. 
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1 Introduction 

The paper concerns problems of long-time convergence to an equilibrium distribution in 
a coupled system. We have proved the convergence for partial differential equations of 
hyperbolic type in IR™ , n > 2 , in [1, 2, 3, 6]; for harmonic crystals in [4, 5], and for a scalar 
field coupled to a harmonic crystal in [7]. Here we treat a particle coupled to a Klein-Gordon 
or wave vector equation. 

Let us outline our main result and the strategy of the proof. (For the formal definitions 
and statements, see Section 2.) Consider the Hamiltonian system consisting of a real- valued 
vector field (f(x) , x G IR 3 , and a particle with position q G IR 3 . The Hamiltonian functional 
is 

TT ( \ 4- f ( \V¥n(x)\ 2 m 2 n \(p n (x)\ 2 \lT n (x)\ 2 x 

H(tp,q,n,p) = 2^ J { 2 2 2 + <Pn(x)q-Vp n (x)) dx 

n=l R3 

+ \{\p\ 2 + u 2 \q\ 2 ). (1.1) 

Here m n > 0, uo > , (p(x) = ((pi(x), . . . , ipd( x )) £ ^ an d correspondingly for n(x) , 
" • " stands for the scalar product in the Euclidean space IR 3 . Taking formally variational 
derivatives in (1.1), the coupled dynamics becomes 

<p n (x, t) = 7r n (x, t), 7r n (x, t) = (A - m 2 n )(p n (x, t) - q{t) ■ Vp n (x), n = l,...,d, (1.2) 

d „ 

q(t) = pit), pit) = -u 2 q(t) -Y, Vn{x, t)Vp n {x) dx, t G IR. (1.3) 

71=1 J 

To state our main results we formulate some assumptions on a constant uj > and a coupled 
function p(x) = (pi(x), . . . , Pd(x)) ■ 

Al. The matrix (u 2 — ml)I — K is positive definite, where m* = if m = (mi, . . . , ma) = 
and m* = min{m n : m n ^ 0, n — 1, . . . , d} if m ^ , K = (Kij)^j =1 is the 3x3 matrix 
with the matrix elements , 

K _v 1 / ^^|p n (A;)| 2 



1 (27r) 3 J k 2 + — m 2 

A2. The function p(x) is a vector real-valued smooth function, p(—x) = p(x) , p(x) = 
for \x\ > R p . 

A3. For all n=l,...,d and k G IR 3 \ {0} , 

p n (k)= je ikx p n (x)dx^0. 

This assumption can be weakened (see Remark 7.6). 

We study the Cauchy problem for the system (1.2)— (1.3) with initial conditions 

<p(x,0)=<p°(x), 7r(x,0) = n°(x), q(0) = q°, p(0)=p°. (1.4) 

Let us write Y = ((p° (x) , q° , tt° (x) , p°) , Y(t) = ((p(x,t),q(t),n(x,t),p(t)) . Then the system 
(1.2)-(1.4) writes as 

Y(t)=F(Y(t)), teWL; Y(0)=Y o . (1.5) 
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We assume that the initial date Y is a random element of a real functional space 8 con- 
sisting of states with a finite local energy, see Definition 2.2 below. The distribution of Y 
is a probability measure /io of mean zero satisfying conditions S1-S3. In particular, the 
correlation functions of the measure /xq is translation- invariant. For a given t e IR, we 
denote by fi t the probability measure on S defining the distribution of the solution Y(t) 
to the problem (1.5). 

Our main result gives the weak convergence of the measures /i t to a limit measure /Xqo , 

Pt^Hoo, t — > oo. (1.6) 

The measure fj,^ is a translation-invariant Gaussian measure on S . Similar results hold as 
t — > — oo because the dynamics is time-reversible. Moreover, in Section 6.3 we prove the 
convergence of the correlation functions of the measures fj, t to a limit as t — > oo . 
We prove the convergence (1.6) by using the strategy of [1, 2] in two steps: 

I. The family of measures fi t , t > , is weakly compact in an appropriate Frechet space. 

II. The characteristic functionals converge to a Gaussian functional, 

fa{Z) = |exp(i(y,Z)) A t t (dy)^exp{-^Q 0O (Z,Z)}, t -> oo, (1.7) 

where Z is an arbitrary element of the dual space and Qoo is the quadratic form defined 
by (2.17), (•, •) stands for the inner product in a real Hilbert space L 2 (IR 3 ) <g) IR^ with 
different N — 1,2, 

Let us explain the main idea of the proof. At first we derive the decay of the order 
(1 + |t|)~ 3 / 2 (and the exponential decay in the case when m = 0) for the local energy 
of the solution Y(t) to (1.5) assuming that the initial date Y has a compact support 
(see Proposition 4.1). Then we apply the integral representation (5.4) of Y(t) and prove 
a uniform bound (5.1) for the mean local energy with respect to the measure pt , t > . 
Finally property I follows from the Prokhorov Theorem. 

To prove of II we derive the asymptotic behavior of the solution Y(t) (see Proposition 
6.1), namely, 

(yW^-EWM/J,^), *-oo, (1-8) 

71=1 

where W n (t) is a solving operator to the Cauchy problem (2.9) for the free wave or Klein 
- Gordon equation, the functions ip^ are expressed by Z (see formula (2.13)). Then we 
apply the results of [1, 2], where the weak convergence of the statistical solutions is proved 
for free wave and Klein-Gordon equations. 

2 Main results 
2.1 Notation 

We assume that the initial data Y are given by an element of the real phase space £ defined 
below. 
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Definition 2.1 Let H = iJ/ oc (IR 3 )©i^ oc (IR 3 ) be the Frechet space of pairs = ((p(x), n(x)) 
with IR d -valued functions (p(x) and n(x) , which is endowed with the local energy seminorms 

H\\l= J (Mx)\ 2 + \V<p(x)\ 2 + \ir(x)\ 2 )dx<oo, Vi?>0. (2.1) 

\x\<R 

Definition 2.2 Let £ = Ti © IR 3 © IR 3 be the Frechet space of vectors Y = (<j>(x),q,p) , 
with the local energy seminorms 

\\Y\\l R =U\\ 2 R + \q\ 2 + \p\\ Vi?>0. (2.2) 

Now we formulate the following condition on p (cf. condition Al). 
Al'. The matrix u 2 I — K is positive definite, where K = (K 0ii j) 3 j =1 is the 3x3 matrix 
with the matrix elements Kq^j , 

kikj\p n (k) \ 2 



„ , (2^./ ~P 



ml 



dk. 



Note that if m — the conditions Al and Al' coincide. If m^O, condition Al' is weaker 
than Al. 

Proposition 2.3 Let conditions Al' and A2 hold. Then (i) for every Yq e £ , the Cauchy 
problem (1.5) has a unique solution Y(t) G C(JR, £) . 

(ii) For every t G IR, the operator U(t) : Y i— > F(t) is continuous on £ . Moreover, for 
every R> R p , T > , 

sup||f/(t)F ||^<C(T)||r ||^ +T . (2.3) 

\t\<T 

Proposition 2.3 is proved in Section 3. 

Let us choose a function ({x) G C °°(IR 3 ) with C(0) + 0. Denote by iJf oc (IR 3 ) , 
s G IR, the local Sobolev spaces of IR d -valued functions </?, i.e., the Frechet spaces of dis- 
tributions ip G -D'(IR 3 ) with the finite seminorms ||v||s,fl := \\^ s {C{ x I R) 1 ^) ||l 2 (r 3 )> where 



A s ip := F k ^ x ((k) s ijj(k)) , (k) := yj \k\ 2 + 1 , and ip is the Fourier transform of a tempered 
distribution ip . 

Definition 2.4 For s G IR denote W = ^"(IR 3 ) © Hf oc (Wi 3 ) , £ s = W © IR 3 © IR 3 . 

Using standard techniques of pseudodifferential operators and Sobolev's Theorem (see, 
e.g., [8]), it is possible to prove that £° = £ C £~ £ for every e > , and the embedding is 
compact. 



2.2 Random solution. Convergence to equilibrium 

Let (Q, E, P) be a probability space with expectation E and B{£) denote the Borel a- 
algebra in £ . We assume that Y = Y (cu,x) in (1.5) is a measurable random function with 
values in (£, B{£)) . In other words, (ou,x) h- >Y (ou,x) is a measurable map 
with respect to the (completed) a -algebra S x B(JR 3 ) and i5(IR 2d+6 ) . Then Y(t) = U(t)Y 
is also a measurable random function with values in (£,B(£)) owing to Proposition 2.3. 
We denote by fio{dYo) a Borel probability measure in £ giving the distribution of Yq . 
Without loss of generality, we assume (f2, S,P) = (£, £>(£), po) and Y (cu,x) = uj(x) for 
Po(duj) x dx -almost all (co,x) G £ x IR 3 . 



3 



Definition 2.5 /i t is a Borel probability measure in £ which gives the distribution of Y(t) : 

Ht(B) = fio(U(-t)B), VBeB(S), teM. 

Our main objective is to prove the weak convergence of the measures Lit in the Frechet 
spaces £~ £ for each e > , 

ih — ? Voo as t -> oo, (2.4) 
where //oo is a limit measure on £ . This means the convergence 
j f{Y)ii t {dY) - j f(Y)^(dY) as t ^ oo 

for any bounded continuous functional /(V) on £~ e . 

Definition 2.6 TTie correlation functions of measure // t are defined by 

Q t (x,y) = E(Y(x,t)®Y(y,t)), for almost all x, y G IR , 

if the expectations in the right hand side are finite. 

We set V = V © IR 3 © H 3 , V : = [C^°(IR 3 ) © IR d ] 2 , and 

(y,Z) := {(j ) ^)+q-u + p-v 

for F = (<f>,q,p) & £ , and Z = (tp,u,v) G £> . For a probability measure // on £ denote by 
fx the characteristic functional (Fourier transform) 

fi(Z) = J exp(i(Y,Z))Li(dY), ZeV. 

A measure \i is called Gaussian (with zero expectation) if its characteristic functional has 
the form 

ft(Z) = exp{~Q(Z,Z)}, ZeV, 

where Q is a real nonnegative quadratic form in P. A measure \i is called translation- 
invariant if 

n(T h B) = li(B), Vfi G B(£), h e IR 3 , 
where T h Y(x) —Y(x — h) . 

2.3 Main theorem 

We assume that the initial measure //q has the following properties SO S3. 

50 //o has zero expectation value, EY (x) = J Y (x) /j,o(dY ) = , x G IR 3 . 

51 iiQ has finite mean energy density, i.e., 

E(\ip°(x)\ 2 + |Vy?°(:r)| 2 + \tc°(x)\ 2 + \q°\ 2 + \p°\ 2 ) < oo. (2.5) 
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Denote by u Q := P/i , where P : (0°, g°, p°) G £ -> 0° G . 
S2 The correlation functions of the measure z/ are translation invariant, i.e., for n,n' G 
J ={!,..., d}, x,yeJR 3 , 



EU{x) ® := C( s - f)> £(Vn(aO ® ^(j/)) := < nn ,(z - y), 

E(n Q n (x) ® <^,(y)) := ^(x - y), ^(s) ® T&fc/)) := - y). 

Now we formulate the mixing condition for the measure vq . 



(2.6) 



Let O(r) be the set of all pairs of open convex subsets A,Bc IR 3 at the distance not 
less than r, d(A, B) > r, and let cr(A) be the a -algebra in Ti generated by the linear 
functionals i— > (0, -0) , where -0 G £>o with supp ip C A. Define the Ibragimov mixing 
coefficient of a probability measure u on H by the rule (cf [9, Def. 17.2.2]) 

/ \ \MAnB)-u (A)u (B)\ 
<f(r) = sup sup — — . (2.7) 

(AB)eO(r) AG (7(^,5 G (7(B) ^ 

z/ (5) > 

Definition 2.7 W^e say i/iai i/ie measure v satisfies the strong uniform Ibragimov mixing 
condition if ip(r) — > as r — > oo . 

S3 The measure v satisfies the strong uniform Ibragimov mixing condition, and 

+oo 

/ r 2 ip 1/2 (r)dr < oo. (2.8) 
o 



Consider the following Cauchy problem 



<p(x,t) = (A-ml)<p(x,t), teTR, 

<p(x,t)\t^o = (p°(x), (p(x,t)\ t=0 = ir (x), x G IR 3 , 



(2.9) 



where m n > 0, ip(x,t) G IR 1 . The following lemma is proved in [1, p. 7], [2, p. 1225]. 

Lemma 2.8 (i) For any 0° = (</?°,7r°) G Tii = iJ/ oc (IR 3 ) © Lf oc (IR 3 ) ? tnere exzsfo a unique 
solution 4>{t) = (<f(x,t),if(x,t)) G C(JR,Hi) to the Cauchy problem (2.9). 
(ii) For any t G IR, the operator W n (t) : 0° i— > 0(t) zs continuous in Hi . 

We define the operator on the space = [Hi] d by the rule 

w(tM, . . . , 0°) = (^i(t)0?, . . . , (2-io) 

Let £ n (x) be the fundamental solution of the operator —A + m 2 n , n G d. For almost all 
x, y G IR 3 , introduce the matrix-valued function Q^^x^y) = (qoo,nn'( x ~ J/))^„, =1 , where 



x ]_ / 9o,ran' + £n * %,nn' %,nn' %,nri \ 

qoo,nn> = (<lZ,nn) ij=Q = Xnn> ~ \ ^ 2 00 ' 

V %,nn> ~ %,nn> %,nn> + ( — ^ + m n)%,nn' J 

where Xnn' = 1 if m n = m n i , and Xnn> = otherwise, the functions g^ nn , > J =0,1, are 
defined in (2.6); and * stands for the convolution of distributions. 
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Remark 2.9 According to [9, Lemma 17.2.3], the derivatives d a qQ nn , are bounded by mix- 
ing coefficient: Va G 2^ with \a\ < 2 — i — j (including a = ), i,j — 0, 1 , 

I^Cn'WI<^ 1/2 (kl), V^GlR 3 , n ,n'6d 

Therefore, <9 a gg n „, G L P (IR 3 ) , p > 1 (see [1, p.16]). Hence, q% nn , G L^IR 3 ) if m n ^ by 
the bound (2.11). If m n = , (2.8) implies the existence of the convolution £ n * q^ nn , . 

Denote by Q^ (i>, 4>) a real quadratic form on S = [S^IR 3 ) ® IR d ] 2 defined by 

d 

Q» 00 ^,^) = (Q» 00 (x,y),^x)®^y)) = E ( goo , nn ,(a; - y), Vn(z) <E> Vv(y)>- (2-12) 

n,n'=l 

The following result can be proved by an easy adaptation of the proof of Theorem B of [1, 2], 
where the result is proved in the case d = 1 . 

Theorem 2.10 Let conditions S0-S3 hold. Then 

(i) the measures v t = W(t)*v weakly converge as t — > oo on £/ie space 7i~ e for each e > . 
(nj T/ie Hmii measure is a translation-invariant Gaussian measure on Ti . 
(Hi) The characteristic functional of is of the form 



Let Z = (ip, u, v) G V , i.e., ip = (V>i, . . . , ipd) G T>$ , (u, t>) G IR 3 x IR 3 . Denote 

d 

i) z = {i>l , . . . , ), := ^fc(^) - E M^) + a fe (x) • u + /3 fc (x) • u, fc G J. (2.13) 

n=l 

Here a k = (a\, a\, a\) , ft = (ft, ft, ft) , 

a£ = <*£(*) = -E/M,( S )^(- S )( VrP fe(x) ^ S , (2.14) 

d +°° / n \ 

# = fi(x) = " E / Mr(s)W£(- S ) ^ v J d a , (2.15) 



M*0 :=E / (^(s)K,n, V>n> 7Z m := I , (2.16) 

i=i £ V WifJn J 

the matrix M(s) = (A/i r (s)) 3 r=1 is defined in (4.17), the operator W^(t) is adjoint to the 
operator W n (t) : 

(<P,W> n (t)^) = (W n (t)<P,^), ip G [S^IR 3 )] 2 , 4> G 7ii, telR. 

Denote by Qoo(Z, Z) a real quadratic form in V of the form 

Q oc (Z,Z) = Q^ z ,^ z ), (2.17) 

where ip z is defined in (2.13). Our main result is the following theorem. 
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Theorem 2.11 Let conditions A1-A3 and SO -S3 hold. Then 

(i) the convergence in (2.4) holds for any e > . 

(ii) The limit measure //oo is a Gaussian equilibrium measure on £ . 
(Hi) The limit characteristic functional has the form 

fi^Z) = expi — Q^Z, Z)}, ZeV. 

(iv) The measure ^ is invariant, i.e., £/(t)*/ioo = /-too, t G IR . 

Remark 2.12 Instead of the strong uniform Ibragimov mixing condition, it suffices to as- 
sume the uniform Rosenblatt mixing condition [13] together with a higher degree ( > 2 ) in 
the bound (2.5), i.e., to assume that there exists a 5 , 5 > , such that 

E (\<p°(x)\ 2+s + |Vy?°(:r)| 2+<5 + \7f°(x)\ 2+s + |g°| 2 + |p°| 2 ) < oo. 

r+oo 

In this case, the condition (2.8) needs the following modification: / r a p (r)dr < oo , 

Jo 

where p = min(<5/(2 + 5), 1/2) , a(r) is the Rosenblatt mixing coefficient defined as in (2.7) 
but without uq(B) in the denominator. 



3 Existence of solutions, a priori estimates 

In this section we prove Proposition 2.3 by the similar arguments as in [10, Lemma 6.3]. 

Let us represent the solution Y(t) as the pair of the functions (Y°(t), Y 1 ^)) , where 
y°(f) = (¥>(*),?(*)), Y\t) = (n(t),p(t)). 

Denote by H S (1R 3 ) the Sobolev space of IR d -valued functions. Let E be the Hilbert 
space of pairs Y = (Y°, Y l ) , where Y° = (</?, q) G H^WL 3 ) © IR 3 , Y 1 = (n,p) G #°(IR 3 ) © 

d 

IR 3 , with the finite norm := ^ (|| Vy?„|| 2 + m 2 n \\(p n \\ 2 + ||7r„|| 2 ) + \q\ 2 + \p\ 2 . Here || • || 

n=l 

stands for the norm in if°(IR 3 ) . Now we prove the auxilary lemma. 
Lemma 3.1 Let conditions Al' and A2 hold. Then 

(i) for every Y G E , the Cauchy problem (1.5) has a unique solution Y(t) G C(IR, E) . 

(ii) For every t G IR, the operator U(t) : Y i— > Y(t) is continuous on E . 
(Hi) The energy is conserved and finite, 

H(Y{t)) = H(Y ) for t G IR. (3.1) 

Proof. Step (i) In the case when p = the existence and uniqueness of solution Y(t) G 
C(IR, E) to the problem (1.5) is proved by Fourier transform. Therefore the problem (1.5) 
for Y(t) G C(IR, E) is equivalent to 

t 

Y(t) = e Aot Y + J e Ao{t - s) BY(s) ds, (3.2) 
o 
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where 



A = 



I 
-Ho 



H Y° = ((-A + m 2 )^, ...,(- A + m 2 )y? d , u; 2 q), 



B 



(Y°,Y 1 ) = (0,RY°), RY :=-(q-V Pl ,...,q-Vp d ,ir U n {x)V p n {x) dx) 

n=l J 



(3.3) 



for Y° = ((/?!, ...,ip d ,q). Note that ||e- 4() *>o||.E < C||^o||b ; and the second term in (3.2) is 
estimated by 

t 

sup || / e A){t - s) BY(s)ds\\ E < CT sup \\Y(s)\\ E . 

\t\<T J \s\<T 

This bound and the contraction mapping principle imply the existence and uniqueness of 
the local solution Y(t) G C([—e,e],E) with an e > 0. 

Step (ii) To prove (3.1) let us assume that 0° = (y?°,7r°) G C 3 (IR 3 ) x C 2 (IR 3 ) and 
<j>° n (x) = for |x| > R , JiG J. Then tp n (x,t) G C 2 (IR 3 x El) and 

ip n (x,t) =0 for |x| > \t\ + max{_R 0) Rp\ 

by the integral representation (3.2) and condition A2. Therefore, for such initial data, the 
equality (3.1) can be proved by integration by parts. Hence the energy conservation (3.1) 
follows from the continuity of U{t) and the fact that [C 3 (IR 3 )] d © IR 3 © [C 2 (IR 3 )] d © IR 3 is 
dense in E . 

Step (Hi) Note that 



>7r) 3 V 



i||Vy? n (:r)|| 2 + im 2 ||v? n (:r)|| 2 + {<p n (x),q- V p n {x)) 



2(2tt) 



k 2 + m 2 n (p n {k) + 



1 



=ig • kp n (k) 



in: 



q ■ k 



III,, 



^n(k)f). 



Hence the Hamiltonian functional is nonnegative, since 



H{Y) = ^E{K|| 2 + 7 ^|| V /fe 2 + m 2 ^(A;)+ 1 iq ■ kp n (k)\\*} 

l Z7r J Jk z + m l n 



n=l 



+l(iU 2 \q\ 2 -q-K q) + hp\ 2 >0 



(3.4) 



by condition Al'. Moreover, by (3.1) and (3.4), for \t\ < e , we obtain that 

\\Y(t)\\ 2 E <CH(Y(t))=CH(Y ). 
On the other hand, for Y = (<p°, n°, q°,p°) , 



^l|vr°|| 2 ) 4 



1 



^)<E(l|V^|| 2 + ^mM| 2 + i 
Hence, we obtain the a priori estimate 

||H0IU<Cill^olU for tern. 



2\|„0|2 



qT + t;\p 



1 2 



(3.5) 



(3.6) 



8 



Therefore, properties (i)-(iii) of Lemma 3.1 for arbitrary t G IR follow from the bound (3.6). 
■ 

We return to the proof of Proposition 2.3. Let us choose R> R p with R p from condition 
A2. Then by the integral representation (3.2) the solution Y(t) for \x\ < R depends only 
on the initial data Y (x) with \x\ < R + \t \ . Thus the continuity of U(t) in £ follows from 
the continuity in E . 

For every R > we define the local energy seminorms by 

imiW=E / (\V Vn (x)f + ml\ Vn (x)\ 2 + \n n (x)\ 2 )dx+\ q \ 2 + \p\ 2 (3.7) 

n=1 \x\<R 

for Y = (ip, it, q,p) . By the estimate (3.6), we obtain the following local energy estimates: 

\\U(t)Y \\ 2 E(R) < C||y |lW|), (3.8) 
for R> R p and t G IR. Hence, the bound (2.3) follows from (3.8). ■ 



4 Decay of local energy 

Proposition 4.1 Let conditions A1-A3 hold and let Y G E be such that 

(p°(x) = tt°(x) = for \x\ > R u (4.1) 

with some R± > . Then there exists a constant C = C(R, Ri) > such that the following 
bound holds for every R > , 

\\Y{t)\\s,R<Ce m {t)\\Y G \\ £ , Rl , t>0, (4.2) 

where 

( _ j e~ 5 ^ with a 5 > 0, if m = , 

£m[t) ~ \ (l + |t|)" 3 / 2 , z/m^0. {A - 6) 

In the case when m = {mi, . . . , 777^) = Proposition 4.1 is extension of Proposition 7.1 
from [10], where a similar result is established in the case when d = 1 and p(x) = p r {\x\) . 
In the case when m^O we apply the methods of [11]. 

To prove Proposition 4.1 we solve the Cauchy problem (1.5) applying the Fourier - Laplace 
transform, 



-oo 



Y(X)= J e- xt Y(t)dt, Re A > 0. 



Then the system (1.2)-(1.4) becomes 



-lfin(x) + \<fn(x, A) = TT n (x, A), U = 1, . . . , d, (4.4) 

-ir° n (x) + \Tt n (x, A) = (A-m 2 n )<f n (x,X)-q(X)-V Pn (x), (4.5) 

-q° + \q(X) = p(X), (4.6) 

- p o + Ap(A) = -uj 2 a(X) -Y,I <Pn(v, A)V Pn (y) dy. (4.7) 

n=l J 
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From (4.4)-(4.7) we obtain 

<p n (x,\) = (-A + ml + XY'i^ + ^-i-A + ml + XY^Pnix)-^), 
(X 2 + u 2 )q(X) = H(X)q(X)+K(X,Y ), (4.8) 

where 

^(A, Y ) := - E /(-A + ml + A 2 ) -1 (<(y) + A^fc/)) V p n {y) dy + (p° + Xq°) , (4.9) 

n=l J 

and H(X) is the 3x3 matrix with the matrix elements Hij(X) , 

d 

Hij(X) = E ^iPn(y)(-A + ml + X 2 )- 1 V jPn (y)dy. (4.10) 

n=l J 

Therefore, equation (4.8) rewrites as 

q(X) = [(X 2 + uJ 2 )I-H(\)Y 1 n(X,Y )=^(\)n(X,Y ), (4.11) 

where by /v (A) we denote a 3 x 3 matrix of the form 

A?(A) = IT 1 (A), £>(A) := (A 2 + cj 2 )/ - #(A) for Re A > 0. (4.12) 



4.1 Time decay for and 

In this subsection we prove the exponential decay for q(t) and pit) in the case when m — . 
The case m ^ is considered in Appendix. 

Theorem 4.2 Lei m = 0, conditions A1-A3 and (4-1) hold. Then there exists a 5 > 
such that the following bound holds 

\q(t)\ + \p(t)\<Ce- St \\Y \\ £ , Rl . (4.13) 

To prove Theorem 4.2 we first investigate the properties of the matrix D(X) . 
Denote (Dp := {X E C : Re A > (3) for /3 G IR. 

Lemma 4.3 Let m = and conditions A1-A3 /ioW. TTien (%) D(X) admits an analytic 
continuation to (D ; (it) for every (3 > 3A^ > such that v ■ D(X)v > C\v | 2 |A| 2 for 
X E (D-^ wjffi |A| > Np and every v E IR 3 . (Hi) There exists a 5 > such that v-D(X)v ^ 
for X E (D _a and /or every r> 7^ . 

Lemma 4.3 is a modification of Lemma 7.2 of [10], where the result is proved in the case 
when d — 1 and p(rr) = p r (M) • 

Proof, (i) We rewrite the entries of the matrix H(X) as 

Hij{X) = E / / V iPn (y) ^ ly _ Zl W jPn (z)dydz. (4.14) 

n=l J J Z \ 



10 



Hence, Hij(X) is denned and has an analytic continuation to (D . Therefore property (i) 
follows. 

(ii) From (4.14) it follows that 

H ij (X) ^ as |A| -> oo with A G (D-p (4.15) 

what implies (ii). 

(iii) At first note that the matrix D(X) is positive definite if ImA = 0. Indeed, from 
condition Al it follows that, for any v G R 3 \ {0} , 

v ■ D(\)v = (X'+^W-it^J^Mktfdk 
> ">l 2 -E^/^l«*>r*>°. 



Secondly, for y G IR 1 \ {0} , we find 



1 f hk.lpMl 2 „ _ + f° r 4 g tJ (r 



Hij(iy + °) - E (27r )3 J k 2 -(y- t0y dk ~ J (r-y + i0)(r + y - iO) ^ 

R 3 



where 

d 



^■(0:=EtAs / Mj\Pn{rO)\ 2 dS 9 . 
™ =1 [Z7r) \e\=i 



) N 

A2. Applying the Plemelj formula (see [12]) and condition A3 we obtain 



Note that gij G C°°([0, +oo)) , max (1 + r) N \gij(r)\ < oo for every N > by condition 

re[0,oo) 



v -1mH(iy + 0)v = -*£j2-pn« I (v ■ 0) 2 \p n (\y\e)\ 2 dS e ^ (4.16) 

1 "=i [Z7l) \e\=i 

for any vector v G IR 3 \ {0} . Hence v • Im D(iy + 0)v — — v • Im H(iy + 0)v ^ 0. Lemma 4.3 
is proved. ■ 

Denote by M(t) an inverse Laplace transformation of N(X) , 

ioo— & 

tf(t) = — f e A 'A/"(A) dX fort>0. (4.17) 

2,711 J 

—ioo— 5 

Lemma 4.4 Let m = and conditions A1-A3 hold. Then for j — 0,1, ... , 

\M {3 \t)\<Ce- & \ t>l. (4.18) 

Proof. By Lemma 4.3 the bound on Af(t) follows. To prove the bound for Af(t) we consider 
XM(X) and prove the following bound: 



v ■ (XAf(X))'v < '. ' , for A G C_«. (4.19) 



C| 




2 


1 + 




2 
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This implies that max \te St J\f(t)\ < oo, and the bound (4.18) for J\f(t) follows. To prove 

te[o,+oo) 

(4.19) it suffices to establish that |A/^(A)| < C(l + |A|)- 3 . From (4.14) it follows that 

H ljW^ Q as \M -»• 00 with Ae(D_ a . 
Therefore, Lemma 4.3 and (4.12) imply that, for i,j = 1,2,3, 

This implies (4.19). The bound (4.18) with j > 2 is proved similarly. ■ 

Proof of Theorem 4.2. Denote by ip^t) = (p\(x,t) the solution of the Cauchy problem 
(2.9) with the initial data 0° = (y?°,vr°) . Then (4.9) and (4.11) imply that for t > 

q(t)=I(t,4> ) + R(t,q°,p ), (4.20) 

where 0° = (0°, . . . , 0°) and 

/(t,0°) := -Y, t f{ip\{t-s)M{s)Vp n )ds 

n=l J Q 
d t 

= E fWs)y^n(t-s),Pn)ds, (4.21) 
n=l|5 

R{t,q°,p°) := M(t) P °+Af(t)q . (4.22) 
First, by the bound (4.18) we have 

\R(t, q ,p )\<Ce- 5t (\q \ + \p \). (4.23) 

Secondly, since m — , 

W n (£)0° =0 for a; G -Br and * > R + R x . (4.24) 

due to the condition (4.1) and the strong Huygens principle. Hence, from (4.18) and (4.21) 
we obtain (see (3.6)) that, for t > R 1 + R p , 

|/(t,0°)| < jtc(p n )f\Af(s)\\\Vti(t-s)\\ LHBRp) ds 

n=l J Q 

< ECi / e-«'|K|| Jll d S <C7 e - a *||^ || Jll . (4.25) 



^ £ — -Rl — 



Relations (4.20), (4.23) and (4.25) imply the bound (4.13) for q(t) . The bound for p(t) is 
proved similarly. ■ 
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4.2 Time decay of field components 

Now we finish the proof of Proposition 4.1. From equations (1.2) we have 

<j> n (x, t) = W n {t)<f n - J W n {t - s) ( ) • ds, n G J. (4.26) 

Let m n = . Then the bound (4.24), the condition (4.1) and Theorem 4.2 imply 

t 

Un(t)\\R <C J e- S9 \\Y \\ e , Rl ds < Ce- St \\Y \\ £ , Rl 

t — Rp — R 

for t > R + max{i? p ,i?!} . Let m n ^ 0. Then instead of (4.24) we apply the following 
well-known bound: 

ii^n(^°ii*<c'(i+r 3/ i0°ik, t>0. 

Hence Theorem 7.1 and the representation (4.26) yield 

\\Mt)\\R<c(i + t)- 3 ' 2 \\Y Q \\ £tRl . 

The bound (4.2) is proved. ■ 

5 Compactness of measures fi t 

Lemma 5.1 Let conditions A1-A3 and S0-S2 hold. Then 

sup E\\U(t)Y Q \\l R < C(R) < oo, V,R>0. (5.1) 



Proof Let us write U (t) = e Aot (see (3.3)). At first, note that 

sup E\\U (t)Y \\ 2 £ R < C(R), \/R>0. (5.2) 

t>o 

Indeed (see (2.1) and (2.2)), 

\\U (t)Y \\l R = \\W(t)<f>°\\l + \go(t)\ 2 + |g (t)l 2 , 
where the operator W(t) is defined in (2.10) and q (t) is a solution to the Cauchy problem 
q {t)+u 2 q {t) = 0, teJR, {qo{t),qo{t))\t=o = {q°,P°)- 

Hence, \q Q (t)\ + \qo(t)\ < C(\q°\ + \p°\) . From [1, Proposition 3.2] and [2, Proposition 3.1] it 
follows that 



d 

I />0 n=\ 



supEWW^Wl = sup 2 E\\W n {t)<tP n f R < C(R), Vi? > 0. (5.3) 
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It implies the bound (5.2). Further, we represent the solution to (1.5) as follows 



U(t)Y = U (t)Y + ju(t- s)BU (s)Y ds, (5.4) 
o 

where the operator B is defined in (3.3). Hence, (4.2) and (5.2) yield 

t 

E\\U(t)Y \\l R < E\\U (t)Y \\l R + Ej\\U(t-s)BU (s)Y \\l R ds 

o 

t 

< C(R) + J e 2 m (t - s)E\\U (s)Y \\l Rp ds < C^R) < oo. ■ 
o 

6 Convergence of characteristic functionals 

6.1 Asymptotic behavior of Y(t) 

Proposition 6.1 Let conditions A1-A3 and S0-S3 hold. Then 

(i) the following bounds hold, 

E\q t (t)-j2( W k(M^l)\ 2 < Ce m (t), (6.1) 
k=i 

E\ Pi {t)-j^(W k {t)^l(3i)\ 2 < Ce m (t), t>l, (6.2) 
k=i 

where e m (t) = e~ 2St , if m — , and e m (t) = (1 + t)^ 1 otherwise, the functions a k , (3\ are 
defined in (2.14) and (2.15), i = 1,2,3. 

(ii) Let ip e [C °°(IR 3 )] 2 with supp ip C Br . Then, for n G d and t > 1 , 



E 



(Mt),1>) - (w n {t)4>l^) + £ (w k {t)<t>le kn ) < ce m (t), 



k=i 



(6.3) 



where the functions 9 kn (x) , k,n e d , are defined in (2.16) with ip n = ip ■ 
Proof, (i) At first, relations (4.20)-(4.22), the bounds (4.23) and (7.11) yield 



d 3 I 

E\qi{t) + E E / (^fc(* " s),Af ir (s)V rPk ) ds 

k =lr=l J 

with e m (t) from (4.3). Secondly, 
E 



< Ce 2 m {t) 



(6.4) 



(fl(t - S), Mr(s)V r Pfc) ^ 

+oo 

J Mr(si) dsi J M ir (s 2 )E ({<pl(t - si), V r p fc )((/4(t - s 2 ),V r p fc )) ds 2 . 
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For any t, s±, S2 G IR , we have 

E((<p 1 k (t-s 1 ),V r p k )(<p 1 k (t-s 2 ),VrPk)) < CsnpE\( V l(r),V rPk }\' 

r6R 



< C iS u P £;||^(t)||| 2(b } <C 2 <oo 
by the bound (5.3). Hence, applying the bounds (4.18) and (7.11) we obtain that 



E 



+00 



(<Pfc(* - s),J\f ir (s)V r Pk) ds 



< Ce m (t). 



Therefore, (6.1) follows from (6.4), (6.5) and (2.14), since 



^(t-s),Mr(s)V r p fc ) = (W k {t)<t>lW' k {-S 



Uir{s)V r Pk 




The bound (6.2) can be proved by similar way. 

(ii) Let i) G [C£°(IR 3 )] 2 with supp^ C B R . From (4.26) it follows that 

3 ' 

(0 n (t), V) = {W n (t)(f>° n , (^(*) K in, ^) ds, 7^ n = 

i=l n 







Note that 

for s > FL p + i? if m n = 0, 
0((l + s)- 3 / 2 ) ifm n ^0. 

Then from (6.1) and (6.7) it follows that 

2 



{w n ( s )n in ,ij) = 



E 



J " s) - Y.(W k {t - s)<f>° k ,ai)) (W n (s)n in ,^) ds 

n k=l 



< Ce m (t). 



Denote 



/„(*) := E\ 



k (t- s)4>° k ,a i k ){W n {s)n in ^) ds 



k=l 



If m n = , /„(*) = for t > + R by (6.7). If m n ^ , 

<Ce m (t). 

It follows from the following estimate: for r G IR , 

d +°° 

S|(^(r)0° fe ,4>r = ^ 



r=l n 



W k (T - s)4>l 



(6.5) 



(6.6) 



(6.7) 



(6.8) 



(6.9) 







ds 



< C < 00, 



by (5.3) and (4.18) and (7.11). The relation (6.6) and the bounds (6.8) and (6.9) imply the 
bound (6.3). ■ 

Corollary 6.2 Let Z = (i/j, u, v) G V = V x IR 3 x IR 3 . Then 

(Y(t),Z) = (W(t)<j>°,i; z )+r(t), 

where (Y(t), Z) = ((p(t),ip) + q(t) ■ u + p(t) ■ v , Y(t) — (<f>(t) , q(t) , p(t)) is a solution to the 
Cauchy problem (1.5), the function ip z is defined in (2.13) and E\r(t)\ 2 < Ce m {t) . 
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6.2 The end of the proof of Theorem 2.11 

Proposition 6.3 Let all assumptions of Theorem 2.11 be fulfilled. Then, for Z G V , 

Eexp{t(Y(t), Z)} - exp { — Q^Z, Z) J , t - oo. 
Proof. By triangle inequality we have 

Ee i(Y(t),z } _ exp{ _I Qoo(Z) z)] < |£( e i<m*> _ e w)^ z >) 



+ 



Se *W*)^>_exp{-i(2oo(Z,Z)} 



(6.10) 



The first term in the RHS of (6.10) is estimated by 



•,{Y{t),Z) _ e i(W{t)^^ z ) 
1/2 



< E 



(Y(t),Z)-(W(t)<f> ,^ 



< E\r{t)\ < (E|r(t)| 2 ) <C4 /2 (t)^0 as t -> oo, 



(6.11) 



by Corollary 6.2. It remains to prove the convergence of E exp{i(W(t)4>°, ip }} = ^(ip ) 
to a limit as £ — > oo . In [1, 2] we proved the convergence of D t (ip) to a limit for ip G Do • 
However, generally, ^ z ^ "Do • Now we introduce a space if m such that ip z G i? m and the 
characteristic functionals i>t{ip) , t G IR, are equicontinuous in if m . 
Denote w n (fc) = (\k\ 2 + m 2 ) 1 / 2 , n G d . 

Definition 6.4 H m is the space of the pairs tp = (ip°, ip 1 ) of IR d -valued functions ip° = 
W,...,^) and V 1 = K,...,^), such that G L 2 (IR 3 ) , < G ^(IR 3 ) , 

n G d , with the finite norm 



E 

n=l 



'.0 1 1 2 



^wMir+iKir+iiv^n 2 ) 



The formulas (2.13)-(2.16) and condition A2 imply that the functions ip z satisfying the 
bound sup fceR 3(l + | | ) 7V (A:) | < oo for every N > , and ip z G H m . Note that if all 
m n ^ , H m = L 2 (IR 3 ) © ff^IR 3 ) . 

Lemma 6.5 (i) The quadratic form Q%(ip,ip) = J | (0°, tp) \ 2 u t (d(f) ) , t G IR, are equicon- 
tinuous in H rn . 

(ii) The characteristic functionals i>t(ip) , t G IR ; are equicontinuous in H m . 



Proof, (i) It suffices to prove the uniform bounds 



sup|Qr(^)l<C||VHL *peH m . 
ten 



(6.12) 



Notethat Q^,V)= E (<lo,nn>(x ~ y), W' n {t)i> n {x) <g> W^(t)^ n ,(y)} . Then by Remark 2.9 

n,n'=l 

we have ^ 

sup |<£(^)| < Csup £ K(t)1> n \\h < CU\\ 2 m , 
ten ten n= \ 
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since in the Fourier transform we have 

F^ k {w n m = ( s ™" n } k)t ) ( f^l ), for i> = 

[ nWrj V u n \k) sm u n {k)t cosuj n {k)t ) \ ^(k) ) ' r , r ; 

(ii) By the Cauchy-Schwartz inequality and (6.12) we obtain that 

I Wi) - MM = I / (^ Vl> - e i <* p -*'>)i/ t (d^ )| < / |£<^-^> - %(#°) 
< / |(0 O ^i-^)k t (# O )<^/|(0 O ,^i-^)| 2 ^(^ ) 

= v/Qr^i - ^2, - < ciivi - Mm- ■ 

Since -?/> z G H m , Proposition 3.3 of [1] (or Proposition 3.2 of [2]) and Lemma 6.5, (ii) 
yield 

v t {i, z ) - exp{-iQ^(^,^)} as t - oo, 
where is defined by (2.12). This completes the proof of Theorem 2.11. ■ 

6.3 Convergence of correlation functions 

Proposition 6.6 Let all assumptions of Theorem 2.11 be fulfilled. Then, for Z X ,Z 2 G V , 

E((Y(t),Z 1 )(Y(t),Z 2 ))^Q O0 (Z 1 ,Z 2 ), t^oo. (6.13) 

Proof. It is enough to prove the convergence of E\(Y(t), Z)\ 2 to a limit as t -> oo . From 
Corollary 6.2 it follows that, for Z G X>, 

E|(F(t),z)i 2 = EK^w/,^)! 2 + o(i) = er(^ z ^ z ) + o(i), t -> oo, 

where ^ z is defined in (2.13) and ^ G . Therefore, Proposition 6.2 of [1] (or Lemma 

4.4 of [2]) and Lemma 6.5, (i) imply that lim^^ Qt(ip z , ip z ) = Q^(ip z , ip z ) . Formula 
(2.17) implies (6.13). ■ 

Acknowledgement 

T.D. thanks Prof. A. Komech for helpful discussions. 

7 Appendix: Time decay of q(t) , p(t) : Case m^O 

Here we prove the following result, which is a modification of Lemma 17.1 of [11]. 

Theorem 7.1 Let m ^ , the conditions A1-A3 hold, and Y G E with (j)°(x) = for 
\x\ > Ri . Then q(t) and p(t) are continuous and the following bound holds 

\q(t)\ + \p(t)\<C(l + tyV\ t>0. 
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Ar first, we rewrite the system (4.4)-(4.7) in the form 

(A - X)Y(X) = -Y , Re A > 0, (7.1) 

where A = Ao + B and the operators Ao , B are defined in (3.3). Hence, the solution Y 
is given by 

Y(\) = -(A- \)- l Y , Re A > 0, 
if the resolvent R(X) — (A — A)" 1 exists for Re A > . 

Proposition 7.2 The operator-valued function R(X) : E —> E is analytic for Re A > . 

Proof. It suffices to prove that the operator (A — A) : E — > E has a bounded inverse 
operator for Re A > . 

Let us prove that Ker (A — A) = for Re A > . Indeed, let Y(X) G E is a solution 
(7.1) with Y = . The function Y(t) = Y(X)e xt e C(JR, E) is the solution to the equation 
Y = AY . Then the Hamiltonian H(Y(t)) = e 2Xt H(Y(X)) grows exponentially by (3.4). 
This grow contradicts (3.1). Hence, (3.4) implies that n (x, X) = 0, p(X) = 0. Equations 
(4.4) and (4.6) imply that <p(x, A) = and q(X) = because A ^ . 

Remember that A = Ao + B (see (3.3)), where the operator B is finite-dimensional and 
the operator A 1 is bounded in E . We rewrite 

A - X = (Ao - X)(I + (Ao - X)- l B), 

where (^4 — A) -1 !? is a compact operator. Since Ker (I + (Ao — A) -1 £?) = , the operator 
/ + (^4 — X)~ 1 B is invertible by the Fredholm theory. ■ 

Denote by gx, n the fundamental solution of the operator —A + m\ + A 2 , 

e -Kn\y\ 

g\, n (y) = . ■ ■ , n e d, (7.2) 

47r|y| 

where n 2 n = m 2 + A 2 , Ren n > for Re A > . From (4.4)-(4.7) it follows that 



M(X) 



9(A) \ _ ( q° 
P(A) ) ~ \ P°(X) 



where p°(A) = p° - E 0/A.n * (A^(A) + n° n (X)), V Pn ) , 

n=l 

M < A ) : =(^/-V) A/)' (7 ' 3) 

and the matrix i?(A) is defined in (4.10). The entries of H(X) are of the form 

HaW = g / V ? p„(y)(, A , n * V jPn )(y) dy = £ JL | dfc. (7.4) 

The following result is proved in [11, p. 351]. 
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Lemma 7.3 (i) The operator — A+m 2 +A 2 is invertible in L 2 (IR 3 ) for ReA > and its 
fundamental solution (7.2) decays exponentially as \y\ — > oo . 

(ii) For every fixed y ^ , the Green function g\ t n(y) admits an analytic continuation (in 

variable A ) to the Riemann surface of the algebraic function \J\ 2 + m 2 n with the branching 
points ±im n if m n ^ . 

Lemma 7.3 and notation (7.3) imply that M(A) admits an analytic continuation from the 
domain Re A > on the Riemann surface with the branching points ±im n with m n ^ , 
n G d. Moreover, the matrix M _1 (A) exists for large A. It follows from (7.3) since 
H(X) -> as ReA -> oo by (7.4). 

Corollary 7.4 (i) The matrix M(A) is invertible for ReA > 0, and 

( 1(A) ) = Ar ' (A) ( /(X) )' ReX > °- (7 ' 5) 

(ii) The matrix M _1 (A) admits meromorphic continuation from the domain ReA > to 
the Riemann surface with the branching points ±im n with m n ^ , n G d . 

Now we investigate the limit values of M _1 (A) at the imaginary axis A = ix , x G IR. 
Without loss of generality, we assume that d = 2 and < m x < m 2 . The other cases can 
be considered similarly. The limit matrix 

M(fe + 0)=( w2/ _^ fa + ()) "J ) , I6 R, (7.6) 

exists, and its entries are continuous functions of x G IR, smooth for \x\ < m\ , m\ < \x\ < 
m 2 , and \x\ > m 2 . 

Lemma 7.5 The limit matrix M{ix + 0) is invertible for x G IR . 

Proof, (i) Let \x\ < mi . Then the matrix (uj 2 — x 2 )I — H(ix + 0) is positive definite. 
Indeed, for every v G IR 3 \ {0} , by the condition Al with m* = m\ , 

f(k-v) 2 \p n (k)\ 2 dk 



v ■ {{uj 2 - x 2 )I - H(ix + 0))v = {uj 2 - x 2 )\v\ 2 - J2{2n)- 3 [ 

n=l J 

> (uo 2 - m\)\v\ 2 - v ■ Kv > 0. 
(ii) Let mi < |x| < m 2 . Then v ■ lmH(ix + 0)v ^ for every v G IR 3 \ {0} . Indeed, 

f . . , f kikj\pi{k)\ 2 dk 

For e > , consider the function 



%(ix + e) 



J k 2 + m 2 — (x — ie) 2 ' 1 



k 2 + m\ — (x — ie) 2 



Denote D e {k) = k 2 + m\ — (x — ie) 2 . For \x\ > mi, D (k) = if \k\ = \Jx 2 — m\. 
We fix a small 5 > and introduce a cutoff function ( G C^°(1R 3 ) such that ((k) > 0, 
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((k) = 1 when \D (k)\ < 5 and ((k) = when \D (k)\ > 25. Note that hnhij(ix + 0) = 
Im hfj(ix + 0) , where 

hUix - 



Denote a(fc) = \Jk 2 + m\. Assume that x > . Since 



1 1 1 

+ 



D e (k) 2a(k)(a(k) -x + ie) 2a(k)(a(k) + x - is) ' 
Im hfj(ix + 0) = Im h s _{ix + 0) , where 

(ia . + £):=/" C(Jb) J^iM dfc . 

We rewrite /ii (ix + e) as 

^(^ + 6)= f-^-da, g(a)= [ ((k) ^[^^ dS. 
~ K J J a + ie ' yv ; 7 sv ; 2a(A;)|Va(A;)| 

a(k)— x=a 

Hence lmh s _(ix + 0) = —ng(0) by the Plemelj formula. Finally, for x > mi > , 

Im h tJ (ix + 0) = -ir f ^y |2 dS. 

|fc| = ^/x 2 — 

Hence, applying condition A3 we obtain that, for mi < \x\ < m 2 , 

v ■ Im H(ix + 0)v = -sign(x)7r(27r)- 3 f dS ^ Q 

\k\=^J x 2 -m\ 

(hi) Let \x\ > m 2 . In this case we find 

2 (27r) -3 j (V ■ kf\Pn{k)? 

n=l 



vlmH(tx + 0)v = -sign(x)vr £ (2tt) 3 / [V ] f nl jl dS ^ 0. 

|fc|=-yx 2 -m2 



In general case c? = 1,2,..., we enumerate mi, ... , m^ in the increasing order, < mi < 
m 2 < . . . < m d . If m fe 7^ m fc+1 and m fc < \x\ < m^+i , 



v!mH(ix + 0)v = -sign(x)7r^(27r)- 3 | ^'^"^^ dg^O (7.7) 

l fe l=-\A 2-m ™ 

by condition A3. For \x\ > m d , the formula (7.7) holds with k = d. ■ 

Remark 7.6 We use condition A3 only in the estimate (7.7). Hence, instead of condition 
A3 it suffices to assume that for any v G IR 3 \ {0} and x > 

£ / (l ,.^|p„(l*|9)| 2 ^)| |th 
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Corollary 7.7 The matrix M 1 (ix + 0) is smooth w.r.t. x G IR outside the points x = 
±imi with mi ^ . 

Lemma 7.8 (%) T/ze matrix M _1 (o;+0) admits the following Puiseux expansion in a neigh- 
borhood of ±imi (mi ^ ): there exists an e± > snc/z £/iai 

oo 

M^tix + O) = J2 c k( x Tm) k/2 , \xTm\<£±, xeJR. (7.8) 

k=0 

(ii) There exists a matrix Ro and a matrix-valued function Ri(x) such that 

M~ l (ix + 0) = -#0 + Ri(x), \x\ > maxm„ + 1, x G IR, (7.9) 

x n 

where \d*Ri{x)\ < Ck/\x\ 2 for \x\ > max n m„ + 1, i e IR, k — 0, 1, 

Proof, (i) Formula (7.8) follows from (7.2) and (7.4). 

(ii) Let / G L 2 (IR 3 ) with supp / C B R . Then (see formula (16.7) of [11]) 

R fe [-A + m 2 + + 0) 2 ]-7Hl W < |x| > m n + 1, 

for every R > . Therefore (see (4.10)) we obtain that \d^Hij(ix + 0)| < Cfc/|x| for 
|x| > max n m n + 1 . Then formula (7.6) implies (7.9). ■ 

Proof of Theorem 7.1. Applying (7.5) we obtain that 

(i§ ) )=h7 e " ,M ^ ix+d Af^+^)) dx - (7 ' io) 

v ' — oo v ' 

Without loss of generality, let us assume that < m\ < m 2 < . . . < m^ . We split the Fourier 
integral (7.10) into d+1 terms by using the partition of unity (o(x) + . . .+Q(x) = 1 , x G IR : 

= e-\C (x) + ... + Ux))M~\ i x + 0)( q ^)dx 

— oo ^ ' 

= I (t) + ... + I d (t), 

where the functions (k(x) G C°°(IR) are supported by 

supp Co C G IR : \x\ > md + 1 or \x\ < mi/2}, 

supp d C {x G IR : mi/3 < |x| < (mi + m 2 )/2}, 

supp ( 2 C {a; G IR : m 2 — 2(m 2 — m±)/3 < |x| < (m 2 + m 3 )/2}, . . . , 

supp Cd C G IR : m d - 2(m d - m d _i)/3 < \x\ < m d + 2}. 

Then (i) the function J (^) G C[0, +00) decays faster than any power of t due to Lemma 
7.8, and (ii) the functions I k (t) G C°°(IR) , k = l,...,d, decay like (1 + |t|)" 3/2 by virtue 
to (7.8). Theorem 7.1 is proved. ■ 

Corollary 7.9 Let m/0 and conditions A1-A3 hold. Then 

W {l \t)\<C{l + \t\y*/\ j = 0,1, (7.11) 

where Hit) is defined in (4-17) and (4-12). This bound can be proved by similar way as 
Theorem 7.1. 
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